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The zero-temperature g-state Potts model partition function for a lat- 
tice strip of fixed width L y and arbitrary length L x has the form P(G, q) = 

J2j={ X c G,j(^G,j) Lx , and is equivalent to the chromatic polynomial for this 
graph. We present exact zero-temperature partition functions for strips of sev- 
eral lattices with (FBC y , PBC X ), i.e., cyclic, boundary conditions. In partic- 
ular, the chromatic polynomial of a family of generalized dodecahedra graphs 
is calculated. The coefficient cqj of degree d in q is c^ = U2d( y ^-), where 
U n (x) is the Chebyshev polynomial of the second kind. We also present the 
chromatic polynomial for the strip of the square lattice with (PBC y , PBC X ), 
i.e., toroidal, boundary conditions and width L y = 4 with the property that 
each set of four vertical vertices forms a tetrahedron. A number of interesting 
and novel features of the continuous accumulation set of the chromatic zeros, 
B are found. 
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I. INTRODUCTION 



The g-state Potts antiferromagnet (AF) |],|2| exhibits nonzero ground state entropy, So > 
(without frustration) for sufficiently large q on a given lattice A or, more generally, on a 
graph G. This is equivalent to a ground state degeneracy per site W > 1, since So = ks In W. 
There is a close connection with graph theory here, since the zero-temperature partition 
function of the above-mentioned g-state Potts antiferromagnet on a graph G satisfies 

Z(G,q,T = 0) PAF = P(G,q) , (1.1) 

where P(G, q) is the chromatic polynomial expressing the number of ways of coloring the 
vertices of the graph G with q colors such that no two adjacent vertices have the same color 
(for reviews, see 0- ||). The minimum number of colors necessary for such a coloring of G 
is called the chromatic number, x(G). Thus 

W({G},q)= lim P(G, q) 1/n , (1.2) 

where n is the number of vertices of G and {G} = lim^oo G. Where no confusion will 
result, we shall sometimes write G rather than {G} for the infinite- length limit of a given 
type of strip graph. At certain special points q s (typically q s = 0, 1, ..,%(£?)), one has the 
noncommutativity of limits 

lim lim P{G,q) 1/n + lim lim P{G,q) 1/n , (1.3) 

q— >q s n— +00 n— +00 q—*q s 

and hence it is necessary to specify the order of the limits in the definition of W({G}, q s ) ||. 
Denoting W qn and W nq as the functions defined by the different order of limits on the left 
and right-hand sides of ( |1.3| ), we take W = W qn here; this has the advantage of removing 
certain isolated discontinuities that are present in W nq . 

Using the expression for P(G, q), one can generalize q from Z + to C. The zeros of P(G, q) 
in the complex q plane are called chromatic zeros; a subset of these may form an accumulation 
set in the n — > oo limit, denoted B, which is the continuous locus of points where M^({G}, q) 
is nonanalytic. [] The maximal region in the complex q plane to which one can analytically 
continue the function W({G},q) from physical values where there is nonzero ground state 
entropy is denoted Ri. The maximal value of q where B intersects the (positive) real axis 
is labeled q c ({G}). This point is important since VF({(j}, q) is a real analytic function from 
large values of q down to q c ({G}). 



For some families of graphs B may be null, and W may also be nonanalytic at certain discrete points. 
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(e) (f) 
FIG. 1. Illustrative strip graphs of the following lattices with length L x = m = 4 and (FBC y , PBC X ) 
boundary condition (cyclic): (a) O m lattice with L y — 4 (b) Q m lattice with L y — 3 (c) S m lattice with 
L y = 3 (d) V m lattice with L y — 3 (e) Xk lT n lattice with k = 3 (f) Yfc, m lattice with fc = 3. 

In the present work we report exact solutions for chromatic polynomials P(G, q) for strips 
of lattices with arbitrarily great length L x vertices and finite width L y . Consider L x sets of 
disjoint tree graphs, T Ly with vertex set v(T Ly ) = {1,2, L y }. Denote the edges joining 
each adjacent pair of T Ly and T' L graphs as e = {viVj\i,j = l,...,L y }, where v^v'j are 
vertices of T Ly and T' L . Impose periodic boundary conditions in the longitudinal direction, 
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denoted x. The lattices in which we are interested here are (a) L y = 4 and e = {11, 32, 44}; 
(b) L y = 3 and e = {11,21,22,33}; (c) L y = 3 and e = {11,21,32,33}; (d) L y = 3 and 
e = {11,21,33}; (e) L y = k and e = {11, kk}; (f) L y = k and e = {11, fcl, fefc}. These 
families of graphs will be denoted as m , Q m , S m ,V m , Xk jJn , Yfc >m in the present work, where 
m = L x . Illustrative examples for these lattices with m = 4 are displayed in Fig. [I]. The 
total number of vertices is L x L y for all these strips. We have given the chromatic polynomial 
for the L y = 3 triangular lattice with e = {11, 21, 22, 32, 33} and cyclic boundary condition 
before [0], and Q m , S m , and V m here are the corresponding lattice strips without one of 
the diagonal edges, without the middle horizontal edge, and without both of these edges in 
each subgraph, respectively. Xk m and Yk^ m come from the homeomorphic expansions on the 
vertical edges of the cyclic L y = 2 strips of the square and triangular lattices. 

A generic form for chromatic polynomials for recursively defined families of graphs, of 
which strip graphs G s are special cases, is [§] 

P((G s ) m ,q) = cg s M{\g s MT , (1-4) 

3=1 

where ca s ,j(q) and the Nq s ,x terms \G s ,j(q) depend on the type of strip graph G s but are 
independent of m. 

The coefficients CG s j{q) that enter into the expressions for the chromatic polynomial ( |1.4| ) 
for the cyclic strip of the families of lattice considered here are 

c (d) = ^(^), (1.5) 

where U n (x) is the Chebyshev polynomial of the second kind, defined by 

[-] 

U n (x)^(-iy( n ~ j )(2xr^ (1.6) 



as we found for the cyclic strips of square and triangular lattices 0. Here in eq. (|1.6|) the 
notation [|] in the upper limit on the summand means the integral part of |. The first few 
of these coefficients are 



c 



(°> = 1, c®=q-l, (1.7) 



.(2) 



3g + l 



and 
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= q 3 -5q 2 + Gq-l . (1.9) 

When the coefficient appears in chromatic polynomials, we shall simply write it as 
unity. 

For a given type of strip graph G s , we denote the sum of the coefficients cg s j as 

C(G S ) = c Gs, • (1.10) 
It will also be convenient to define the following polynomial: 

^^^tI-b-^C; 1 ) 91 " 2 "*' (L11) 

where P(C n , q) — (q — l) n + (q — 1)(— l) n is the chromatic polynomial for the circuit (cyclic) 
graph C n with n vertices. Some works on calculations of chromatic polynomials for recursive 



families of graphs include JTT|- [45[] 



II. FAMILY OF O m STRIPS WITH (FBCy,PBC x ) 

In this section we give our solutions for the chromatic polynomials of the 4 x m strips of 
the O m lattice strip with cyclic boundary conditions. The chromatic number for this family is 
x(0) = 3. We calculate the chromatic polynomials by iterated use of the deletion-contraction 



theorem, via a generating function approach |p^ , p| , and a coloring matrix method ||13|| . We 
find N ,x = 14 and 

14 

P(O m ,FBC y ,PBC x ,q) = J2 c oA*o tj ) m , (2.1) 

3=1 

where Aoj's for 1 < j < 13 are roots of a cubic equation, a sixth-degree equation, and a 
quartic equation. Specifically, the Aoj's for 1 < j < 3 are roots of the equation 

£ 3 - (<? 4 - 6g 3 + I6q 2 - 24q + 17)^ 2 

+(q -l)(q- 2)(g 4 - 7q 3 + 20q 2 - 28q + 18)^ - (q - l) 2 (q - 2) 2 = (2.2) 
while the Ao,/s for 4 < j < 9 are roots of the equation 
£ 6 + 2(q 3 - 6q 2 + Uq - 13)^ 5 + (q 2 - 5q + 7)(g 4 - 9q 3 + 30q 2 - 47q + 29)^ 4 
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-(2g 8 - 30g 7 + 202g 6 - 796g 5 + 2001g 4 - 3278g 3 + 3420g 2 - 2084g + 567)£ 3 

+(g - l)(g 9 - 17g 8 + 130g 7 - 586g 6 + 1717g 5 - 3401g 4 + 4578g 3 - 4072g 2 + 2199g - 557)£ 2 

-2(g - l) 2 (g - 2)(g 5 - 9g 4 + 33g 3 - 61g 2 + 57g - 23)£ 

+ (g-l) 4 (g-2) 2 = 

the Aoj's for 10 < j < 13 are roots of the equation 

£ 4 - (q 2 - 6g + 10)£ 3 - (2g 3 - 13g 2 + 30g - 23)£ 2 

-(q - l)(g - 3)(g 2 - 4g + 5)£ + (g - l) 2 = (2.4) 

and finally, Ao,i4 = 1. 

The corresponding coefficients are 

c OJ = 1 for 1 < j < 3 (2.5) 

Co j = c (i) for 4 < j < 9 (2.6) 

c OJ = c ( 2 ) for 10 < j < 13 (2.7) 

and 

co,i4 = c (3) . (2.8) 

Here the coefficients have the form of the Chebyshev polynomial of the second kind given in 
( |1.6|) -( |1~9[ ). There is also a zero eigenvalue Ao,i5 = with coefficient Co, 15 = q(q — l)(g 2 — 
3g+ 1), which does not contribute to the chromatic polynomial ( p.l|) . The sum of all of the 
coefficients is equal to F(T 4 ,g) = g(g — l) 3 if we include Co, 15, and is equal to g 2 (g — 1) if 
Co, 15 is not included. 

The locus B separates the g plane into three main regions. The outermost one, region 
Ri, extends to infinite |g| and includes the intervals g > g c and g < on the real g axis. 
Here, 



g c (0, FBC y , PBC X ) ~ 2.638342 



(2.9) 



Region R2 includes the real interval 2 < q < q c , while region R3 includes the real interval 
< q < 2. In regions i?,, 1 < z < 3, the dominant terms are, respectively, the root with 
maximal magnitude of (1) the cubic equation (|2.2j ), (2) the quartic equation (p.4j) , (3) the 
sixth-degree equation (|2~3l ). Our previous calculations for various families of graphs p0|,|7| 
have shown that B can include pairs of extremely small complex-conjugate regions. We have 
not made an exhaustive search for these in the present case. The locus B for this family 
(defined in the limit n — > 00) and, for comparison, chromatic zeros for a typical long finite 
O m strip, with L x = m = 20, are shown in Fig. [| With this value of m, the chromatic zeros 
lie close to the boundary B. The locus B crosses the real g-axis at q = 0, 2 and q c . 



Im(q) 




Re(q) 
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FIG. 2. The locus B for the m — > oo limit of the O m family of cyclic strip graphs (solid curve). For 
comparison, we show chromatic zeros for the cyclic O m strip graph with m = 20, i.e., n = 80 vertices. The 
dashed line is the locus of solutions of the degeneracy equation involving only the roots of the cubic eq. 



Earlier, it was shown that the locus B can have support for Re(q) < in |3(J for square 



lattice strips and in [[/]] for triangular lattice with cyclic or Mobius boundary conditions. 
These families have respective girths equal to 4 and 3, respectively, for L x large enough to 
avoid degenerate cases. Chromatic zeros with negative real parts were also found for several 
other families of graphs |f[7| , 24] , including those that yielded chromatic zeros with arbitrarily 



large magnitude \q\ |p6| , |33|] . Here we find that B(0) also has support for Re(q) < for the 
present m — > oo limit of) the O m family, which has girth 5. Although the chromatic roots 
of regular dodecahedron, denoted as O5 here, does not have any negative real part, where 
the smallest one is about 0.0082126, except for the trivial root q = 0, the chromatic roots 
of O m for m > 6 do include some with negative real parts. The largest magnitude of the 
negative real part increases as m increases, and approaches the limit of B up to q ~ —0.125. 
The density of zeros is observed to be relatively small on the parts of B extending through 
q = and 2, separating R% from R± and R2 near the real g-axis. 

Partial results for the chromatic polynomial of the O m family were reported in ||44j| , where 



this family was denoted as D n . Where they overlap, some of the results in our complete 



calculation differ with those in jPJ E6 



III. FAMILY OF Q M STRIPS WITH (FBC Y ,PBC X ) 

In this section we consider the chromatic polynomials of the 3 x m strips of the Q m 
lattice with (FBC y , PBC X ), i.e. cyclic, boundary conditions. This family of graphs has one 
more diagonal edge in each subgraph than the L y = 3 square lattice with cyclic boundary 
condition p9| , and one less diagonal edge than the L y = 3 triangular lattice with the same 



boundary condition 0. We have Aq A = 10, and 
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P(Q m ,FBC y ,PBC x ,q) = J2 c Q>i( X Qj) m > ( 3 - X ) 

where 

^0,(1,2) = \ (V - 6g 2 + Uq - 13 ± y/(q - 3)(g 2 - 5q + 7)(q* - V + 6q - 5) ) . (3.2) 
The Aqj for 3 < j < 6 are the roots of the quartic equation 
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£ 4 + (3g 2 - 15g + 20)£ 3 + (3g 4 - 30g 3 + 114g 2 - 191g + 117)£ 2 



+ (g - 2)(g 5 - 14g 4 + Teg 3 - 199g 2 + 251g - 121)^ - (q - 2) 5 (g 2 - 5g + 5) = (3.3) 



Aq, (7 ,8) = g - 3 ± ^4 - g (3.4) 
and \q$ = q — 2, Aq^o = —1. The coefficients are 

c QJ = 1 for 1 < j < 2 (3.5) 

c Q j = c (1) for 3 < j < 6 (3.6) 

CQ j = c ( 2 ) for 7 < j < 9 (3.7) 



and 



c Q ,io = c (3) . (3i 



Again only one type of polynomial of each degree in q occurs in the coefficients, and 
the number of A's with coefficients c^ d \ which is denoted as np(L y , d) in 0, are exactly the 
same as those for the square and triangular lattices with L y = 3: np(3,0) = 2,np(3, 1) = 
4,np(3,2) = 3,np(3,3) = 1. Therefore, the sum of the coefficients is 

10 

C(Q, FBCy, PBC X ) = CQ,j = l(q - l ? ■ (3-9) 

3=1 

For the family Q m of lattice strips, we find that q c = 3; at this point all the Aq's have the 
same magnitude 1. The locus B crosses the real g-axis at g = 0, 2 and g c = 3, and separates 
the g plane into five main regions. Region Ri extends to infinite |g| and includes the intervals 
g > 3 and q < on the real q axis. Region R 2 includes the real interval 2 < q < 3, and region 
i?3 includes the real interval < q < 2. There is also a pair of complex-conjugate regions 
i? 4 , i?4 centered at approximately q ~ 2.55 ± 1.30i. We have not made an exhaustive search 



for the tiny regions in this case. The dominant terms in regions Ri and R 2 are Aq^ in (3.2) 
and Aqj in (|3.4j ), while regions R3 and R4, RI are both dominated by the roots of the quartic 
equation ( |3~3"1) with maximal magnitude. The locus B has support for Re(q) > 0, and the 
density of zeros is smaller on the parts of B extending through q = and 2, separating R3 
from Ri and R 2 - 

The locus B and chromatic zeros for the L x = m = 20 cyclic graph of the Q m lattice are 
shown in Fig. 13. The chromatic zeros lie close to the boundary B and indicate its position. 
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FIG. 3. Locus B for the m — > oo limit of the family of cyclic graphs Q m - For comparison, chromatic 
zeros for Q m with m = 20 (i.e., n — 60) are shown. 

IV. FAMILY OF S M STRIPS WITH (FBC Yl PBC x ) 

The 3 x m strips of the S m lattice with (FBC y , PBC X ) = cyclic boundary conditions 
have one less middle horizontal edge in each subgraph than the L y = 3 triangular lattice 
with the same boundary conditions [j/j. For the solutions of the chromatic polynomial, we 
find N S) \ = 6 and 

6 

P(S m , FBCy, PBC X , q) = J2 c sA*sj) m , (4-1) 

3=1 
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where 





) 



(4.2) 





) 



(4.3) 



and A 



■s,5 — 



(q — 2) 2 , Xs t Q = q — 4. The coefficients are 



C S'>(1>2) 



1 



(4.4) 



c s,j 



c (1) for 3 < j < 5 



(4.5) 



and 



CSfi = c 



(4.6) 



These coefficients have the form of Chebyshev polynomials of the second kind, as given 
above. There is also a zero eigenvalue Xsj = with coefficient csj = q(q 2 — 3q + 1) which 
does not contribute to the chromatic polynomial (|4.1| ). The sum of the coefficients is equal 
to P(T 3 , q) = q(q — l) 2 if we include csj, and is equal to q 2 if cgj is not included. 
The locus B crosses the real g-axis at q = 0, 2 and q c , where 



which is the real solution of equation 2q 3 — llq 2 + 24g — 24 = 0. This locus separates the q 
plane into three regions. Region Ri extends to infinite \q\ and includes the intervals q > q c 
and q < on the real q axis. Region R 2 includes the real interval 2 < q < q c , while region 
i?3 includes the real interval < q < 2. The dominant terms in regions R4, 1 < % < 3 are (i) 
A^i in (|4.2D , (ii) A^^, and (iii) Xs,3 i n C N=. 3|) , respectively. Therefore, q c for the m — >• 00 limit 
of this lattice strip S m is given by the degeneracy between IA^jJ and \Xs,q\, which yields the 
value in ( [4.71 ). The locus B and the chromatic zeros for the L x = m = 20 cyclic graph of the 
S m lattice strip are shown in Fig. [|. This locus has support for Re(q) > 0. The locus B(S) 
not only has closed regions R\ to R3, but also has a pair of arcs extending from q ~ 2.9±0.7i 
at the boundary between R 1 and R 3 to q ~ 3.05 ± 0.396i. 



q c (S, FBCy, PBC X ) = 2.702873 



(4.7) 
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FIG. 4. Locus B for the m — > oo limit of the cyclic family SVn. For comparison, chromatic zeros for the 
graph S m with m = 20 (i.e., n = 60) are shown. 



V. FAMILY OF V M STRIPS WITH (FBC Y ,PBC X ) 

We next proceed to the 3xm strips of the V m lattice with (FBC y , PBC X ), which have one 
less horizontal edge in each subgraph than the Q m lattice, or equivalently one less diagonal 
edge in each subgraph than the S m lattice. We find AV,a = 6 as for lattice S m , and 

6 

P(V m , FBC y , PBC X , q) = J2 c vj(Vj) m , (5.1) 

3=1 

where 
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\v,(i,2) = ^ (f ~ Sq + 4 ± y /( q -2)(q-3)(q*-q + 2) ) , (5.2) 
Xvj for j = 3 < j < 5 are roots of the cubic equation 
£ 3 + (2g 2 -9g+ll)£ 2 

+(g - 2)(g 3 - 8g 2 + 20q - 14)£ - (g - l) 2 (g - 2) 2 (q - 3) = (5.3) 
and Ay,6 = g — 3. The corresponding coefficients are 

cy,(i,2) = 1 (5.4) 

c VJ = c « for 3 < j < 5 (5.5) 

and 

cy, 6 = c (2) . (5.6) 

These coefficients are the same as the coefficients for the family of lattice strips S m , and the 
zero eigenvalue Xyj = has the same coefficient cyj = q(q 2 — 3q+l). Thus, the sum of the 
coefficients is equal to F(T 3 , q) = q(q — l) 2 if we include cy,7, and is q 2 if cyj if not included. 

The locus B has support for Re(q) > 0, and consists of (1) an outer closed region which 
crosses the real (/-axis at q = 0, |; (2) an inner closed region which crosses the real g-axis 
at q = 2 and q ~ 2.341; and (3) a finite line segment on the the real g-axis extending from 
q ~ 2.341 to q = 3. Thus, the g plane is separated into three regions and q c = 3. Region Ri 
is the exterior part of the outer closed region, and extends to infinite \q\. Region R3 is the 
interior part of the inner closed region, and region R2 includes the interior part of the outer 
closed region and exterior part of R%. The dominant terms in regions Ri, 1 < i < 3 are (i) 
Ay ( i in ( |5.2| ), (ii) the roots of eq. ( |5T3| ) with largest magnitude, and (iii) Ay,6, respectively. The 



locus B for the m — ► 00 limit of the V m family and the chromatic zeros for the L x = m = 20 
V m graph are shown in Fig. |5|. As was true for the S m lattice, The locus B(V) has both 
closed regions and finite line segment. 
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FIG. 5. Locus B for the m — > oo limit of the cyclic V m family. For comparison, chromatic zeros for the 
cyclic V m graph with m = 20 (i.e., n = 60) are shown. 



VI. FAMILY OF X K ,M STRIPS WITH (FBC Y ,PBC X ) 
The homeomorphic expansion of the L y = 2 square lattice with cyclic boundary condi- 



tion on the horizontal edges has been studied before [27]. In this section we consider the 
homeomorphic expansion of the L y = 2 square lattice with cyclic boundary condition on 
the vertical edges, so that there are k vertices on the transverse direction, i.e. we make the 
insertion of k — 2 degree-2 vertices on all of the vertical edges of the strip. This is equivalent 
to the L y = k square lattice having the same boundary condition without all the internal 
horizontal edges. We find N x ,\ = 6 as for lattice V m , and 
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P(X k , m ,FBC y ,PBC x ,q) = Y t '-sJ'\s i ,jY" , (6-1) 

3=1 



where Xx k ,(i,2) are the roots of the quadratic equation, 

:-i) k 

q 



e+ ( ~^- (V - 2q + 2)(1 - q)™ + q - 2^ 



q 2 



((g - l) (2fc - 3) - (q - 2)(1 - g) (fc ~ 2) - l) = . (6.2) 
The Ax fe ,(3,4) are the roots of another quadratic equation, 

e 2 + ^((2-g)(l-g) (fe - 1) -2)e 

-- 2 ({q - l)^ 1 ) - (g - 2)(1 - q) k - (q - l) 2 ) = (6.3) 



and 



X Xkl5 = ~(q - l)D h {q) , (6-4) 



X Xkfi = D k {q) . (6.5) 



where D k {q) is defined in In eqs. ( |6.2|) and (|6.3| ), the coefficient functions are 

polynomials for any positive integer k. 
The corresponding coefficients are 

cx,(i,2) = 1 (6.6) 
c XJ = c (i) for 3 < j < 5 (6.7) 

and 

cx,e = c {2) . (6.8) 

These coefficients are the same as the coefficients for the lattices S and V. Therefore, the 
sum of the coefficients is equal to q 2 if we do not consider the zero eigenvalue. The chromatic 
number is 2 and 3 when L x is even and odd, respectively, independent of k. 

This lattice is the only one we consider here which does not contain a seam when we 
impose periodic longitudinal boundary conditions with an orientation-reversal (Mobius). 
We have Xx k ,mb,j = Xx k j, 1 < j < 6, and the corresponding coefficients are 
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Cx,mb,(l,2) — 1 



(6.9) 



c X,mb,(3,4) — c 



(6.10) 



c X,mb,5 — 



(6.11) 



and 



CX,mb,6 — — 1 • 



(6.12) 



The changes of the coefficients when we change the longitudinal boundary condition from 
cyclic to Mobius follow the rules given in Theorem 6 of 0. The chromatic number for the 
Mobius strips of lattice X k , m is 2 as (k,L x ) = (o, e),(e,o), and 3 as (k, L x ) = (o, o),(e,e), 
where e and o denote even and odd values for k and L x . 

The fact that the X/s for a Mobius strip must be the same as those for the cyclic strip of 
the same width and lattice type was proved in [[31]]; this also proves, a fortiori, that (i) the 
total number, ]Va, of A/s, and (ii) the continuous nonanalytic locus £>, including the point 
q c , are the same for the cyclic and Mobius strips. 
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FIG. 6. Locus B for the m — ► oo limit of the cyclic X3, m family and chromatic zeros for the cyclic ^3, 
graph with to = 20 (i.e., n = 60). 
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FIG. 7. Locus B for the m — > oo limit of the cyclic Xi. m family and chromatic zeros for the cyclic Xi^ m 
graph with m — 20 (i.e., n = 80). 

The locus B and the chromatic zeros for the L x = m = 20 cyclic graph of the X^ iJn lattice 
with k = 3, 4 are shown in Fig. |] and [7]. Except for k = 2 which is the ladder graph, we find 
the locus B(Xk) for k > 3 has support for Re(q) < 0, and crosses the real g-axis at q = 
and 2. The two main regions are region R\ includes the intervals q > 2 and q < on the real 
g axis, and region R2 includes the real interval < q < 2. The dominant terms in regions 
Ri and R 2 are roots of (|6.2|) and ( |6.3|) , respectively. 

For odd k, there is also a pair of complex-conjugate regions R 3 , R^ on the upper right and 
lower right of q = q c = 2, which are connected to R2 by arcs. The dominant term in these 
regions is Xx odd ,3- The locus B(X dd) around q c has the form of a circular arc which bends 
to the left before reaching a pair of complex-conjugate roots of the discriminant of eq. (|6.3|) . 
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For even k, a general feature is that there are two contiguous pairs of complex-conjugate 
regions R 3 , on the upper right and lower right of q = 2 and _R 4 , R\ on the upper left and 
lower left of q = 2, and both of them touch the real g-axis at q = 2. The dominant terms in 
i? 3 ,i?3 and R 4 ,R* 4 are Xx even ,5 and Ax e „ en ,6, respectively. 

As k increases, more small regions appear on the boundary between regions Ri and R%. 
These small regions surround part of the roots of the discriminant of eq. ( |6.3|) , as shown in 
Fig. H for k = 7 and Fig. |9] for k = 8. In these figures the roots of the discriminant are 
indicated with an asterisk, *. Note that the locus B is compact and approaches the circle 
\q — 1| = 1 as k increases, as is illustrated in Fig. |10| for k = 15. This is similar to what 
was found before for the outer envelope of the locus B for the homeomorphic expansions on 
the horizontal edges of the cyclic L y = 2 strips of the square lattices [p7| . The roots of the 
discriminant of eq. (|6.3|) also approach the same circle as k increases, as is illustrated in Fig. 



11 for k = 50. 
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FIG. 8. Locus £> for the m — > oo limit of the family of cyclic graphs -XV,rn and chromatic zeros for the 
cyclic graph Xj^ m with m — 15 (i.e., n — 105). The roots of the discriminant of eq. (pit) are indicated with 
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FIG. 9. Locus B for the m — > oo limit of the family of cyclic graphs Xg,m and chromatic zeros for the 
cyclic graph X%^ m with m — 15 (i.e., n — 120). The roots of the discriminant of eq. d6.3| ) are indicated with 
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FIG. 10. Locus B for the m — > oo limit of the family of cyclic graphs X\$^ m and chromatic zeros for the 
cyclic graph X\^ iTn with m = 6 (i.e., n = 90). 
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FIG. 11. The roots of the discriminant of eq. (3.3) with k = 50. The solid circle is \q — 1| = 1 



VII. FAMILY OF Y KM STRIPS WITH (FBC Y ,PBC X ) 

In this section we insert k — 2 degree-2 vertices on all the vertical edges of the L y = 2 
triangular lattice with cyclic boundary conditions. Here we find iVy,A = 4, and 



P{Y k>m , FBCy, PBC X , q) = J2 °rA X Y k 



j) ' 



(7.1) 



where 



X Yk>1 = (q-l)D k+1 (q)-(q-2)D k (q) 
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= D,(q)D k (q) + (q-l)(-l)^ , 



(7.2) 



The Ay fci ( 2 ,3) are the roots of the quadratic equation, 
e 2 + td£( 2 (2-g)(l-#- 1 ) + g -4)c 

+^ ((<? - 3)(g - l)^*- 1 ) -(q- 2)(2q - 3)(1 - q)^ + q 2 - 3q + 3 ) = , (7.3) 
where the coefficient functions are polynomials for positive integer k, and 

\ YkA = D k (q) . (7.4) 
The corresponding coefficients are 

cy,i = 1 (7.5) 
cy,( 2 ,3) = c (1) (7.6) 

and 

cy, 4 = c< 2 > . (7.7) 

These coefficients are the same as the coefficients for the square and triangular lattices 
with L y = 2 and the same boundary condition. Therefore, the sum of the coefficients is equal 
to P(T 2 ,q) = q(q — 1) if we do not consider the zero eigenvalue. The chromatic number is 
2 as (k,L x ) = (o, e), and 3 as (k,L x ) = (o, o), (e,e), (e,o), where e and o denote even and 
odd values for k and L x . The locus B for the m — > oo limit of the cyclic family Yfc jm and the 



chromatic zeros for the L x = m = 20 cyclic graph Yk t m with k = 3,4 are shown in Fig. [12 
and 113 
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FIG. 12. Locus B for the m — > oo limit of the family of c^clfc graphs Ys iTn and chromatic zeros for the 
cyclic graph Y3 ;TO with m = 20 (i.e., n = 60). 
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FIG. 13. Locus B for the m — > oo limit of the family cRe^C})c graphs l4, m and chromatic zeros for the 
cyclic graph Y^ m with m = 20 (i.e., n = 80). 

The locus B has support for Re(q) > 0, and crosses the real g-axis at q = 0, 2 for odd 
fc and g = 0,2, g c for even k, where 2 < q c < 3. The two main regions for odd k are region 
which includes the intervals q > 2 and g < on the real q axis, and region i?2, which 
includes the real interval < q < 2. For even k, there is another region R3 which includes 
the real interval 2 < q < q c , together with region f?i including the intervals q > q c and q < 
on the real q axis. The value of this q c decreases as k increases. The dominant terms in 
regions R\ and R 2 are Ay i in (|7.2j ) and the root of eq. ( |7.3|) . For even k, the dominant term 
in R 3 is Ay j4 in (|7.4j ). 

For k = 3, there is also a pair of complex-conjugate regions R3, Rl contiguous to R2, and 
the dominant term in these regions is Ay 3i 3. The B{Yk) around q c again has the form of a 
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circular arc that bends to the left before reaching a pair of roots of the discriminant of eq. 

(OD- 

For k > 4, there are more regions that appear inside R2; in these regions, the dominant 
term is Xy k ,i- These internal regions are connected to an arc, and the free end of the most 
right-hand arc is connected to the boundary between regions Ri and i?2- The most left- 
land arc has two small regions on both of its ended if the small region is not connected 
to the boundary between regions R\ and R2. The free ends of these arcs are roots of the 
discriminant of eq. (|7.3|) , and the small regions they are connected to also surround the 
roots of the discriminant. This is illustrated in Fig. |14] and Fig. [15], where the roots of the 
discriminant are indicated with * . These internal regions approach the circle | q — 1 1 = 1 as 



k increases, as is indicated in Fig. [16] for k = 15. The roots of the discriminant of eq. (|7] 
also approach the same circle as k increases, as is illustrated in Fig. [17] for k = 50. 
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FIG. 14. Locus B for the m — > oo limit of the cyclic fanB^^E}-),,,. and chromatic zeros for the cyclic graph 
^7,m with m — 15 (i.e., n — 105). The roots of the discriminant of eq. ([7.3[) are indicated with *. 
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FIG. 15. Locus B for the m — > oo limit of the cyclic fanB$ (E|) TO and chromatic zeros for the cyclic graph 
Ys t m with m — 15 (i.e., n — 120). The roots of the discriminant of eq. ( [7.3[ ) are indicated with *. 
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FIG. 16. Locus B for the m — > oo limit of the cyclic famB^|^ !m and chromatic zeros for the cyclic graph 
5^15, m with m — 6 (i.e., n — 90). 



29 



1.5 



0.5 - 



- 



-0.5 



-1.5 




-0.5 0.5 1 1.5 



2.5 



FIG. 17. The roots of the discriminant of eq. (7.3) with k = 50. The solid circle is \q — 1\ = 1 



VIII. COMPARATIVE DISCUSSION 



We have found some interesting features of the locus B for these strips with 
(FBCy, PBC X ) = cyclic boundary condition: 

1. For the families of the lattice strips that we have studied so far, as long as the boundary 
conditions are cyclic, i.e., each repeated vertex set forms a tree graph (Tl ) and the 
last vertex set is connected to the first vertex set to have periodic boundary conditions 
in the longitudinal direction, the coefficients CG sJ (g)'s that enter into the expressions 
for the chromatic polynomial (|1.4j) are always in the form of Chebyshev polynomial 
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of the second kind given in eq. ( |1.5|) . This property was found to be true for square 
and triangular lattices with cyclic boundary condition in |J. Note that, while this 
statement applies also to certain classes of nonplanar strip graphs, such as the Mobius 
strips of the square lattice, it does not apply to the nonplanar strip graphs that we 



have analyzed in for these, we showed that the coefficients have different forms 
and are similar to the coefficients of the strips with doubly periodic boundary condition 
(torus) 0,H|5 . 



2. In general, the sum of the coefficients C(G S ) in eq. (|1.10|) is the total dimension of the 
space of coloring configurations and should be equal to P(Ti y ) = q(q — l) L y _1 for the 
lattice strips with (FBC y , PBC X ) boundary conditions. Lattice Q m does satisfies this 
relation, as shown in eq. (|3.9|) . For other lattice strips studied here, the sum of the 
coefficients without taking account the coefficient of the zero eigenvalue are not P(T Ly ), 
but this can be understood easily. Consider two adjacent vertex sets of lattice X k>m 
with k > 3 and apply the coloring matrix from one vertex set to another. Effectively, 
only the two end vertices are connected and should be considered for coloring. Since 
they are separated by other vertices in between, the total number of proper colorings 
of these vertices is q 2 . A similar argument applies to the lattice strip V m in Fig. [I] 
(d), and the lattice strip S m in Fig. [18] (a). For the lattice strip O m , effectively only 
two connected vertices (forming a T 2 ) and a separated vertex in each vertex set are 
connected to another vertex set; therefore the sum of the coefficients is q 2 {q — 1). With 
the homeomorphic expansions on the vertical edges of the cyclic L y = 2 strips of the 



triangular lattice, Yk,m does not satisfy the condition C given in ||43 |. However, it is 



equivalent to the homeomorphic expansions on the diagonal edges of the same lattice 
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strip, as shown in Fig. [T^ (b), so the sum of the coefficients is still q(q — 1). 





3. The loci B's for the infinite-length limits of the strips studied here with global circuits in 
the longitudinal direction cross the real q-axis at q = 0, 2 as before. The feature of the 
cyclic strips that there is a correlation between the coefficient cg s j of the respective 
dominant A<? s j's in regions that include intervals of the real axis discussed in |7j is 
found to be true again. The ca a js of the dominant A^j's in (1) region Ri containing 
the intervals q > q c and q < 0, (2) region containing the intervals < q < 2, and (3) 
region containing the interval 2 < q < q s , where q s may not be q c , are = l,c^\ 
and S- 2 \ respectively. 

4. In previous calculations of B for the zero-temperature Potts antiferromagnet on infinite- 
length limits of strip graphs with global circuits (i.e., PBC X ) P,pO|. 
the respective loci did not involve arc endpoints. Here we find such arcs with endpoints 
such that one end is connected to the boundary of regions for the infinite-length limits 
of the families S, V, and with k > 7, and arcs connecting two regions for lattices V, 
X odd and Y k with k > 4. In Table | we list various properties of these lattices presented 
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in this work. 



5. Although some loci B calculated for the infinite- length limits of families of graphs with 
global circuits involved disconnected pieces (e.g., p6jp7| ), here we find a series of such 



disconnected components, with a remarkably intricate structure, as is evident in Figs. 



14, 15, and 16 



TABLE I. Properties of P, W, and B for strip graphs G s of lattices studied here. The properties apply 
for a given strip of type G s of size L y x L x ; some apply for arbitrary L x , such as N\, while others apply for 
the infinite-length limit, such as the properties of the locus B. The boundary conditions is (FBC y , PBC X ) 
= cyclic. The column denoted eqs. describes the numbers and degrees of the algebraic equations giving 
the Ag sJ -; for example, {2(1), 2(2), 1(4)} indicates that there are 2 linear equations, 2 quadratic equations 
and one quartic equation. The column denoted BCR lists the points at which B crosses the real q axis; the 
largest of these is q c for the given family G s . The notation "int;qi;g c " refers to cases where B contains a 
real interval, there is a crossing at qi, and the right-hand cndpoint of the interval is q c . q c for lattice Y^ m 
with even k depends on the value of k, and is not shown here. Column labeled "SN" refers to whether B 
has support for negative Re(q), indicated as yes (y) or no (n). 



G s 


Ly 


N x 


eqs. 


BCR 


SN 





4 


14 


{1(1)1(3)1(4)1(6)} 


0, 2, 2.638 


y 


Q 


3 


10 


{2(1)2(2)1(4)} 


0, 2, 3 


n 


s 


3 


6 


{2(1)2(2)} 


0, 2, 2.703 


n 


V 


3 


6 


{1(1)1(2)1(3)} 


0, 2, 2.341, int: 2.5; 3 


n 


X 


k 


6 


{2(1)2(2)} 


0, 2 


y 


Y 


k 


4 


{2(1)1(2)} 


0) 2, Qc(^ei)en) 


n 



IX. CONCLUSIONS 



In this paper we have presented exact calculations of the zero-temperature g-state Potts 
antiferromagnet partition function (equivalently, chromatic polynomial P), for various lattice 
strips of fixed width L y and arbitrarily great length L x with cyclic boundary conditions 
(FBCy, PBC X ). The coefficient cqj of degree d in q is = U^di^)' wnere U n (x) is the 
Chebyshev polynomial of the second kind. We explain the sums of the coefficients for the 
respective families. We also find a number of interesting and novel features of the singular 
locus B. The results herein are thus of interest both from the viewpoint of exact results in 
statistical mechanics, graph theory, and algebraic geometry. 
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X. APPENDIX: FAMILY OF Z M STRIPS WITH (PBC Y ,PBC X ) 



In this section we consider the 4 x m strips of the square lattice with (PBC y , PBC X ) 
boundary conditions (torus), but with two diagonal edges in each basic vertex set , i.e., 
instead of L x set of tree graphs, we have L x set of complete graphP] K4 (tetrahedron) to 
begin with and e = {11, 22, 33, 44}. The lattice with length L x = m will be denoted as Z m , 
and an illustrative example with m = 4 is displayed in Fig. |l£. 
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FIG. 19. Illustrative strip graphs of Z m lattice with L y = 4 with length L x = m = 4 and (PBC y , PBC X ) 
= torus boundary condition. Vertices are indicated with • (points where edges cross without a symbol • are 
not vertices.) 



In ||43|| , the lattice with L x set of for arbitrary L y is called the bracelet strip, and 
the A's with coefficients 1, O(q), and 0(q 2 ) are given. For L y = 2, it is the ladder graph 
with vertex set = K2. For lattice strip with L y = 3, which has vertex set C3 = K3, the 
chromatic polynomial was given in [3^|. Here we give the full set of A's for the L y = 4 strip. 
We find N z \ = 16 and 
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P(Z m ,PBC y ,PBC x ,q) = J2 c zA X zj) m . 

i=i 



(10.1) 



where all A^j's for 1 < j < 16 are polynomials: 



2 K n denotes the complete graph, i.e. the graph with n vertices such that each vertex is connected 
by an edge to all of the other vertices. 
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\ Z1 = q 4 - I0 q 3 + Uq 2 - 84q + 73 (10.2) 

A Zi2 = -(g 3 - 12q 2 + 50g - 73) (10.3) 

\z,3 = -(q-3)(q 2 -5q + 7) (10.4) 

A Zi4 = q 2 - 9g + 21 (10.5) 

Az, 5 = 4 2 -5g + 5 (10.6) 

A z , 6 = g 2 -llg + 31 (10.7) 

A z>7 = 3 2 -5g + 7 (10.8) 

Az, 8 = 9 2 -7g+ll (10.9) 

Az,9 = i-g (io.io) 

Az,io = 2-g (10.11) 

Az,n = 3-g (10.12) 

Az,i2 = 4-g (10.13) 

Az,i3 = 5-g (10.14) 

Az,i4 = 6-g (10.15) 

Az,i5 = 7-g (10.16) 
and Xz,i6 = 1- The corresponding coefficients are 

cz,i = 1 (10.17) 
35 



and 



cz,2 = \c z ^ = q-l (10.18) 



c ZA = c Zi5 = ^(q-l)(q-2) (10.19) 



111.. . 
cz,e = 2°zj = 3 C ^,8 = ^Q{Q ~ 3) (10.20) 



cz,9 = \c z , 13 = \{q - l)(q - 2)(q - 3) (10.21) 
3 6 



cz.io = c ZM = ^c z , 12 = q{q - 2)(q - 4) (10.22) 



cz,i 5 = \c z ,u = \q{q - - 5) (10.23) 
3 6 



c z l6 = g 4 - 10g 3 + 29g 2 - 24g + 1 . (10.24) 



The sum of all of the coefficients is equal to P(C 4 ,g) = q(q — l)D±(q). The chromatic 
number is x{Z~) = 4. 

The locus B and chromatic zeros for the L x = m = 20 cyclic graph of the Z m lattice are 
shown in Fig. The locus B crosses the real (/-axis at q = 0, 2 and q c , where 

q c (Z, PBC y , PBC X ) ~ 3.673593 . (10.25) 
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FIG. 20. Locus B for the to — > oo limit of the family Z m with toroidal boundary conditions and chromatic 
zeros for Z m with to = 20 (i.e., n = 80). 

The locus B has support for Re(q) > 0, and separates the q plane into three regions. The 
outermost one, region R\, extends to infinite \q\ and includes the intervals q > q c and q < 
on the real q axis. Region f?2 includes the real interval 2 < q < q c , while region R% includes 
the real interval < q < 2. In regions Ri, 1 < i < 3, the dominant terms are Xzi, ^z,e, 



and Xz,2, respectively. Thus, the q c {Z) given in (|10.25|) is the degeneracy between \\z,x\ and 
\Xz&\i an d is the rea l solution of q 3 — 9q 2 + 31g — 42. 



The locus B has three regions and crosses the real g-axis at q 
quite similar to the B for the same lattice strip with L y 



0,2, and q c , which is 



3 |32| but has a larger size. Since 
the dominant A's in these three regions have coefficients with order 0(1), O(q), and 0(q 2 ), 
which are called level 0, 1, and 2 in [43j, naively one may infer this to be true for arbitrary 
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Ly. The locus B determined from the eight general A's given in [[|3] for 3 < L y < 10 i 



is 



shown in Fig. |21], where the envelope of the outer boundary of B does increase in size as 
L y increases. However, for L y = 8,9, and 10, an extra small regions R± appears around 
q = 7, 8, and 8, respectively. Therefore, even for such simple structure of A's (in the form of 
polynomial), the locus B may not remain simple for sufficiently large L y . 
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FIG. 21. Loci 23's for the infinite-length limit of the square lattice strip with K^ y as the basic subgraph 
vertex set and (PBC X , PBC y )~ torus boundary conditions. In order of the sizes of the main envelopes from 
small to large, the 23's are for 3 < L y < 10. 



Acknowledgment: I would like to thank Prof. R. Shrock for helpful discussions. 



38 



R. B. Potts, Proc. Camb. Phil. Soc. 48 (1952) 106. 
F. Y. Wu, Rev. Mod. Phys. 54 (1982) 235. 
R. C. Read, J. Combin. Theory 4 (1968) 52. 

R. C. Read and W. T. Tutte, "Chromatic Polynomials", in Selected Topics in Graph Theory, 
3, eds. L. W. Beineke and R. J. Wilson (Academic Press, New York, 1988.). 
N. L. Biggs, Algebraic Graph Theory (2nd ed., Cambridge Univ. Press, Cambridge, 1993). 
R. Shrock and S.-H. Tsai, Phys. Rev. E55 (1997) 5165. 



S.-C. Chang and R. Shrock, cond-mat / 0005236 



S. Beraha, J. Kahane, and N. Weiss, J. Combin. Theory B 27 (1979) 1; ibid., 28 (1980) 52. 
S.-C. Chang and R. Shrock, Physica A, in press (|cond-mat/0005232|) . 
E. H. Lieb, Phys. Rev. 162 (1967) 162. 

N. L. Biggs, R. M. Damerell, and D. A. Sands, J. Combin. Theory B 12 (1972) 123. 
N. L. Biggs and G. H. Meredith, J. Combin. Theory B20 (1976) 5. 
N. L. Biggs, Bull. London Math. Soc. 9 (1977) 54. 

R. C. Read, in Proc. 3rd Caribbean Conf. on Combin. and Computing (1981). 
R. C. Read, in Proc. 5th Caribbean Conf. on Combin. and Computing (1988). 
R. J. Baxter, J. Phys. A 20 (1987) 5241. 

R. C. Read and G. F. Royle, in Graph Theory, Combinatorics, and Applications (Wiley, NY, 
1991), vol. 2, p. 1009. 

R. Shrock and S.-H. Tsai, Phys. Rev. E55 (1997) 6791. 

R. Shrock and S.-H. Tsai, Phys. Rev. E56 (1997) 1342, 2733, 4111. 

R. Shrock and S.-H. Tsai, Phys. Rev. E56 (1997) 3935. 

J. Salas and A. Sokal, J. Stat. Phys. 86 (1997) 551. 

M. Rocek, R. Shrock, and S.-H. Tsai, Physica A252 (1998) 505. 

M. Rocek, R. Shrock, and S.-H. Tsai, Physica A259 (1998) 367. 

R. Shrock and S.-H. Tsai, Physica A259 (1998) 315. 

R. Shrock and S.-H. Tsai, Phys. Rev. E58 (1998) 4332; |cond-mat/9808057| . 

R. Shrock and S.-H. Tsai, J. Phys. A 31 (1998) 9641; Physica A265 (1999) 186. 

R. Shrock and S.-H. Tsai, J. Phys. A Lett. 32 (1999) L195. 

R. Shrock and S.-H. Tsai, J. Phys. A 32 (1999) 5053. 

R. Shrock and S.-H. Tsai, Phys. Rev. E60 (1999) 3512. 

R. Shrock and S.-H. Tsai, Physica A 275 (2000) 429. 

R. Shrock, Phys. Lett. A261 (1999) 57. 

N. L. Biggs and R. Shrock, J. Phys. A (Letts) 32, L489 (1999). 
A. Sokal, |cond-mat/9904146 . 



R. Shrock, in the Proceedings of the 1999 British Combinatorial Conference, BCC99 (July, 
1999), Discrete Math., in press. 
[35] R. Shrock, Physica A281, 221 (2000). 



39 



N. L. Biggs, reports LSE-CDAM-99-03,05. 

S.-C. Chang and R. Shrock, Phys. Rev. E62, 4650 (2000). 

R. Shrock, Physica A283, 388 (2000). 

S.-C. Chang and R. Shrock, Physica A286, 189 (2000). 

S.-C. Chang and R. Shrock, Physica A, in press (|cond-mat/0004161| ). 

S.-C. Chang and R. Shrock, Phys. Rev. E62, 4650 (2000). 

S.-C. Chang and R. Shrock, |cond-mat/0008477| . 

N. L. Biggs, LSE report LSE-CDAM-00-04 (May 2000). 

N. L. Biggs and P. Reinfeld, LSE report LSE-CDAM-00-07 (June 2000), 

[http://www.cdam.lse.ac.uk/Reports/ . 



S.-C. Chang and R. Shrock, Physica A, in press ( cond-mat / 000749 1| ) . 

The differences include the chromatic polynomial, the singular locus, and the question of 
whether there are chromatic zeros with Re{q) < 0. We are informed by Prof. Biggs and P. 
Reinfeld that they now agree with our calculation and conclusions and are revising LSE- 
CDAM-00-07 accordingly. We thank Prof. N. Biggs and P. Reinfeld for correspondence on 
this. 



40 



